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The bands of graphite are extremely sensitive to topological defects which modify the electronic 
structure. In this paper we found non-dispersive flat bands no farther than 10 meV of the Fermi 
energy in slightly twisted bilayer graphene as a signature of a transition from a parabolic dispersion 
of bilayer graphene to the characteristic linear dispersion of graphene. This transition occurs for 
relative rotation angles of layers around 1.5° and is related to a process of layer decoupling. We have 
performed ab-initio calculations to develop a tight binding model with an interaction Hamiltonian 
between layers that includes the tt orbitals of all atoms and takes into account interactions up to 
third nearest-neighbors within a layer. 



Graphene is a sheet of carbon atoms arranged on a honeycomb lattice recently obtained by micro-mechanical 
cleavage of graphite^. This two-dimensional lattice has a singular linear band dispersioiP which makes the charge 
carriers behave as massless Dirac fermions with a speed of Vp s» 10 6 m/s, travelling large distances without interaction. 
These properties change drastically in the presence of a second layer. A Bernal or AB stacked bilayer graphene (BLG) 
is a stack of two carbon sheets in such a way that an atom of one of the layers is in the center of the other layer's 
hexagon. Like single layer graphene, BLG is also a scmimctal but its dispersion relation is quadratic and its charge 
carriers have a non-zero effective mass. 

Few layers of graphene grown epitaxially on different surfaces show a variety of defects including rotations of the 
top layer, for instance in graphene layers prepared by chemical vapor deposition (CVD) 3 . Moire patterns are very 
often obs erv ed in Scanning Tunneling Microscopy (STM) measurements when a relative rotation between top layers 
is presenlPS 

Graphite is also extremely sensitive to topological defects which can modify its electronic structure. Extended 
defects like lattice dislocations lead to the presence of localized states at the Fermi energj03 ag j s t ne case f graphene 
ribbons with zigzag edges^. These localized states can also be found as a result of local defects such as in graphene 
antidot lattices^. They were predicted in superstructures with honeycomb symmetry by N. Shima et alP^, who suggest 
the occurrence of ferromagnetism when electron correlation is turned on. 

In addition one of the mechanisms proposed to explain high- T c superconductivity is associated with the presence 
of extended Van Hove Singularities (VHS) near the Fermi energyMSl. This kind of V HS ar ising from a nearly dis- 
persionless band has been observed in high-T c cuprates by angle-resolved photoemissiorP^Ell. The superconductivity 
in graphene when the fermi level is close to a VHS has been also explored theoreticalljff^. Recently Guohong Li 
et alP^ reported the observation of two symmetric low-energy VHSs in the density of states, measured by scanning 
tunneling spectroscopy, in twisted graphene layers. They showed that the position of these singularities can be tuned 
by controlling the relative angle between layers. 

In this paper we predict the presence of flat bands, close to the fermi level, in slightly twisted BLG for angles around 
1.5°. We show that for rotation angles below 1°, as expected, the parabolic band dispersion of BLG is obtained, while 
for angles above 10° up to the symmetric point of 30°, the dispersion is linear as in single layer graphene at the 
K pointfi^ni. For rotation angles between 1° and 2° our results show a transition to a linea r disp ersion, otherwise 
between 2° and 10°, a linear dispersion prevails with a renormalized velocity of Dirac fermiond 20 * 21 [ 

A rapid flattening of the bands is observed when we start rotating a stacked AB Bilayer, flat bands very close to 
the Fermi level are revealed for angles around 1.5°. This effect has been reported in Refp2J were they observe two 
pronounced peaks in the DOS measured by STM, and by changing the rotation angle between the layers they achieve 
a constant shift of the VHS's toward Fermi level. In this work, based on a band structure analysis, we show the 
existence of a critical angle for which the energy separation between the VHS's reaches a minimum. 

We have developed a tight-binding (TB) model to tackle commensurable unit cells with a large amount of atoms. 
The model considers the interlayer interaction of tt orbitals of all atoms between layers and up to third nearest- 
neighbors interaction within a layer. We fit our TB model to reproduce the band structure obtained from Density 
Functional Theory (DFT) calculations of a stacked AB bilayer. 

For building a commensurate unit cell we have followed a procedure proposed by Kolmogorov and CrespP^l. We 
started our rotations from a stacked AB bilayer. In this stacking there are two non-equivalent sites, A (a) site, where 
an atom lies directly above another atom and, B site where an atom position is just in the center of another layer 
hexagon. We have chosen a B site as our rotation center, however the results are essentially the same if we choose an 
A site, the difference is that with the latter we can obtain a stacked A A superstructure (all atoms in one layer are 
above/below another). 

We do a clockwise commensurable rotation from a vector r = ma\ + na,2 to t\ = na\ + ma2, where a\ = 
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FIG. 1: (Color Online) Energy band dispersion for 6 different structures. Twisted BLG (straight blue line), single layer 
graphene (red dotted) and DFT results (green diamonds). At the Fermi level, around the if -point a linear behavior of the 
bands is observed for all structures. In e and f the Fermi velocity is renormalized. 



(— 1/2, y/3/2)ao and S2 = (1/2, -\/3/2)ao are bilayer lattice vectors; m,n are integers and ao = 2A6A is the lattice 
constant. The unit cell vectors are : t\ = na\ -\-ma2 and fa = —ma\ + (n + m)(Z2 and it contains ./V = 4(n 2 + mn + ra 2 ) 
atoms. 

We have adopted the tight binding model proposed by Reich et alPHl to study the electronic structure of graphene, 
which includes third-nearest neighbor approximation for the non-orthogonal 7r orbitals model. Such parameterization 
properly describes the 7r and 7r* graphene bands and it is able to reproduce first principle results over the entire Brillouin 
zone. The Hamiltonian for twisted stacks is given by H = Hi + H2 + Hi n t, where Hi and H2 are the Hamiltonians for 
each layer following the Reich model and H{ n t = r y t e~( r ~ d ^ 13 ■ is the interaction Hamiltonian between the two layers. 
Here d is the interlayer distance and j t and f3 are parameters. 

By choosing this interaction Hamiltonian we are taking into account the complexity of the unit cell where the 
distances between atoms of different layers are all different. There are nine parameters in this TB model and they 
were obtained fitting the band structure of a stacked AB bilayer to reproduce DFT results. In our calculations the 
interaction between layers is not restricted to the nearest neighbors, we have included interaction between all atoms 
of different layers. If only the nearest neighbor is included in the interlayer interaction, the AA site gets a higher 
significance (a unique site in each cell where one atom lies exactly above another), overestimating the renormalization 
of fermi velocitjEH. 

Ab-initio calculations of the electronic structure of BLG and some twisted configurations (namely 
(2,1);(3,2);(4,3),(5,4)) were done using the SIESTA^ code with atomic basis DZP (double-^ with polarization), us- 
ing the LDA^ approximation for the exchange correlation potentials with an energy cut of E cut = 270R y . The 
Monkhorst Pack grid 10 x 10 x 1 was used as a k-points grid. All structures were fully relaxed using the Conjugate 
Gradient algorithm CG with force tolerance of 0.02eVyA. The interlayer distance used was c = 3.35A We have 
compared the band structures obtained by SIESTA calculations with those obtained with our TB model for the four 
rotated stacks mentioned, we found no difference in the low energy range. In FigjT] (b),(c) and (d) it is shown the 
corresponding band structures for different cells calculated with DFT and with the TB model. 

In figure [T] we show the behavior of twisted bilayer graphene and single layer graphene for angles between 3.9° to 
nearly 30°, in six different configurations. The K and M points have different values for each cell, so for comparison 
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FIG. 2: (Color online) A flattening of the bands is observed with a depletion in vf for the slightly twisted Bilayer (blue line). 
The AB stacked BLG (green dashed line) and single layer graphene (red dotted line) are included for comparison. 

we included the one layer and/or bilayer graphene energy bands in every plot, which are calculated with the same 
unit cell vectors of the twisted bilayer graphene. 

At the Fermi level and around the K point, a linear band dispersion of the twisted cell can be seen in all graphics. 
In Fig(lja), (b), (c) andJIJd), the slope of the dispersion curves is identical to those of a single layer graphene and in 
FigjTJe) andjljf) a decreasing of the band slope is evident around the K point in the rotated BLG. 

For rotation angles above 10° the weak coupling between layers is responsible for this behavior. The Fermi velocity 
vf at the K point of the rotated bilayer is the same as in single layer graphene in agreement witlP^HSD An increase 
in the interlayer distance, obtained by DFT relaxation 2 -, only enhances this effect. It is important to notice that in 
the review published by Pong et alp3 most of the Moire patterns have been observed for angles below 10°. These 
structures do not produce Moire Patterns strong enough to be detected by STM^l 

In figure [2] the differences in the slope of the twisted structure compared with that of the single layer graphene are 
more evident. The stacked AB BLG bands are included to enhance visual comparison. It is also clear that bands 
become less dispersive for lo w angles and they touch the T point closer to the Fermi level. Our results are in agreement 
with those obtained b} * 18 * 21 ! using different approaches. 

For relative rotation angles of 1.6° and 1.3° (Fig[3ja),(b)) the bands are completely flat. The inset in[3jL) reveals 
a near zero Fermi velocity near K ponrpSl. By comparing the insets of Fig|3] (a) and FigJJ^b) a change in the band 
convexity and the presence of asymmetric flat bands very close to the Fermi level can be observed. Otherwise, for an 
angle of 1° (Fig|3jc) the behavior is parabolic as expected. 

Since BLG has a parabolic dispersion, we would not expect a drastic change of the band structure if we twist it 
by a very small angle. By increasing the angle, the layers start a decoupling process and this influences the band 
dispersion. Around 1.5° a transition characterized by flat bands occurs. The energy of these bands changes with the 
angle and so does the relative extension of the flat area. These flat bands are only present for a small range of angles 
and the separation between them, near the Fermi energy, reaches a minimum at 6 C = 1.5°. The relative extension of 
the flat bands is maximized at the same angle. This latter behavior can be associated with localized states near fermi 
level. 

The behavior shown by the flat bands mentioned in the preceding paragraph is also observed for the energy gap 
AE (energy difference between the minimum of the conduction band and the maximum of the valence band) at the T 
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FIG. 3: (Color online) Localized states in the band structure of twisted BLG for angles around 1.5°. Flat bands are not 
present for 1°, instead a parabolic dispersion appears. Insets in a and b are obtained near the K point. They show a change 
in the convexity of the band from T to K point and, in b, a near zero Fermi velocity is revealed. 



point. In Fig[3|b) we have plotted this energy gap as a function of the rotation angle. It is clearly seen in the graphic 
the existence of a critical angle that separates two different behavior of this energy gap. We have fitted the results 
by separating the data above and below the critical angle 9 C = 1.5° and we can observe a linear behavior of AE as a 
function of 9 for angles greater than 9 C and a quadratic behavior, for smaller angles. This transition point indicates 
two qualitatively different electronic behavior, a massive particle for small rotation angles between the graphene layers 
and massless Dirac fermions when the angle exceeds the critical value. A slight shift in the value of the critical angle 
is quite possible provided the descriptive character of the TB results. 

The fact that the behavior of AE at T point is an indication of the presence of flat bands in Twisted BLG is very 
convenient because we only need to calculate one point (r) of the band structure to plot this graph, otherwise it 
would be very time consuming since unit cells for angles below 1° contain more than 15000 atoms. 

In FigQa) we also plotted the quotient between the Fermi velocity in twisted BLG and the corresponding Fermi 
velocity in single layer graphene, v F /v F , as a function of the twisting angle. By fitting the data we have found 
that this velocity is depleted for an angle of ~ 1.5°. For smaller angles a parabolic dispersion it is retained and thus 
regarding the presence of Dirac electrons it is not adequate. They are only present for angles above this critical angle. 

We believe that the arising of these flat bands in the twisted structure is related to a lost of a degree of freedom in 
the system, when layers decouple the electrons are confined to each layer, this transition would cause the appearance 
of extended VHS. 

In conclusion we have studied the band structure of a slightly twisted bilayer graphene by employing ab-initio 
calculations to develop a parameterized tight binding(TB) model. This model has allowed us to deal with large 
unit cells and also to compare our results with those previously reported. Our Hamiltonian includes all interactions 
between atoms on different layers and takes into account interactions up to third nearest-neighbors within a layer. We 
have found the presence of flat bands at the K point around a critical rotation angle of 1.5°. This is the signature of a 
transition from a parabolic to linear dispersion in the twisted structure. A degree of freedom is lost in this transition, 
layers decouple and electrons are confined to each layer for angles above the critical one, and the system behaves 
like separated graphene layers. For this critical angle the Fermi velocity in a twisted BLG is completely depleted. 
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FIG. 4: (Color online) a. Dependence of Fermi velocity in twisted BLG v F and in single layer graphene v F with the rotating 
angle, b. AE at T point versus the rotation angle. 



Therefore, the presence of Dirac fermions in twisted bilayer graphene is valid only for angles larger than the critical 
one. Why this transition occurs precisely at this angle of 1.5° is probably dependent upon our parametrization model, 
but we have found that a critical angle certainly exists. 
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